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ABSTRACT 

The aim of this work is to test to what extent the star cluster population of a galaxy can be 
utilised to constrain or estimate the star formation history, with the Large Magellanic Cloud as 
our testbed. We follow two methods to extract information about the star formation rate from 
star clusters, either using only the most massive clusters (following Maschberger & Kroupa 
|2007| l or using the whole cluster population, albeit this is only possible for a shorter age span. 
We compare these results with the star formation history derived from colour-magnitude dia- 
grams and find good overall agreement for the most recent « 1 Gyr. For later ages, and espe- 
cially during the "cluster age gap", there is a deficiency of star clusters in relation to the star 
formation rate derived from the colour-magnitude diagram. The star formation rates follow- 
ing from the whole cluster population lie a factor of s» 10 lower than the star formation rates 
deduced from the most massive clusters or from the colour-magnitude diagram, suggesting 
that only sa 10% of all stars form in long-lived bound star clusters. 

Key words: galaxies:evolution - galaxies: individual (Large Magellanic Cloud) - galaxies: 
stellar content - galaxies: star clusters 



1 INTRODUCTION 

The understanding of galaxy evolution is a major goal of astro- 
physics. Every large-scale event in the life of a galaxy, as e.g. an in- 
teraction with another galaxy, has its own pattern of star formation. 
Since stars can have long lifetimes, the stellar population preserves 
information of such events, allowing one to re-trace the galaxy's 
evolution from the present stellar content. In this work we focus on 
the star formation history, the progression of the star formation rate 
in time. This study has two main aspects, the comparison of two 
different methods to obtain a star formation history, from colour- 
magnitude diagrams and from the star cluster population using the 
method of Maschberger & Kroupa ( 2007 1. Furthermore we discuss 
the star formation history of the Large Magellanic Cloud, which 
serves as a "guinea pig" for the comparison. 

The common method to obtain a detailed star formation his- 
tory is to observe the stars in a galaxy (or a part of it). From the 
distribution of the stars in a colour-magnitude diagram (CMD) the 
star formation rate at a given time can be derived using modelled 
tracks of stellar evolution. To get a result which is representative 
for the whole galaxy it is necessary to observe a significant frac- 
tion of the stars in the galaxy, distributed over a large area. This 
leads to limitations of this method: since individual stars need to 
be resolved, only nearby galaxies can be examined. Also, a large 
number of stars and a large area demand a big observational effort. 
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Fortunately, the Large Magellanic Cloud has been extensively ob- 
served, so that a set 24 million stars is available, from which [Harris| 
|& Za ritsky (2009) derive the star formation history. 

Another approach to infer a star formation history was pre- 
sented by Maschberger & Kroupa (2007). Here the fact is used 
that practically all stars form in star clusters. The notion of a star 
cluster is here taken in a wider sense, denoting stellar assemblies 
from the smallest size, say a dozen stars, up to classical globular 
clusters, and does not necessarily imply a bound system. Whereas 
low-mass clusters will disperse their stars rather quickly into the 
galactic field, massive clusters have lifetimes comparable to a Hub- 
ble time. Since the time distribution of massive clusters is related 
to the star formation rate at their birth, they can be used to find the 
star formation history of their host galaxy. This approach using the 
most massive star clusters has been investigated by Maschberger & 
|Kroupa| ( |2007| l from a theoretical point of view, applying Monte- 
Carlo models to study how reliable the massive clusters trace the 
star formation history. The value of this method is that, as (mas- 
sive) star clusters are observable up to much further distances than 
individual stars, the star formation history of a wider range and 
number of galaxies can potentially be obtained. 

Additionally the time-sequence star formation occurring in 
star clusters can be determined by simply taking all star clusters 
into account, given that also a large fraction of lower-mass star 
clusters is usually observed. Subramaniam (2004) compared in six 
regions of the Large Magellanic Cloud the cluster formation rates 
and star formation rates (derived from colour-magnitude diagrams), 
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Figure 1. Overview of the observed regions in the Large Magellanic Cloud 
(background image from Bothun & Thompson 1 19881 with astrometry by 
Parker et al.||1998J. The dots mark the centres of the fields observed by 



Harris & Zaritsky ( 2009 1 and for which a star formation history was derived 
using a colour-magnitude diagram. The squares are the boundaries of the 
regions observed by Massey ( 2002 1, in which |de Grijs & Anders] ^2006) 
derived star cluster ages and masses. The open dots are the fields of |Harris| 
& Zaritsky 1 2009 1 which we selected for comparison. 
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Figure 2. Age-Mass diagram of the star clusters in the Large Magellanic 
Cloud(grey circles original data from de Grijs & Anders 2006 1. The cluster 
masses have been corrected for dynamical evolution using eq.[T]with U = 
8 Gyr (black dots, top panel) and = 1 Gyr (black dots bottom panel). 
30Dor, not part of the original data set, is shown as a star. 



and found that in general the number of formed clusters follows the 
"field" star formation rates. Cluster formation rates only have also 
been derived by|Girardi et al.|l|1995),|Pietrzynski & Udalski| ([ 2000| > 



|Hunter etal (2003), and de Grijs & Anders {2006} Our approach is 
slightly different, we use the total mass in clusters per time instead 
of number of clusters per time. With the diminishing brightness of 
star clusters as they age the observed number of older star clusters 
decreases, so that only a shorter fraction of the galaxy's life-time 
can be investigated in this way. 

The Large Magellanic Cloud has been the target of many re- 
search projects and is ideal to compare the CMD and star cluster 
methods. In the literature there are a number of studies available on 
the star formation history of the Large Magellanic Cloud, which are 
based on the CMD approach (e.g. Harris & Zaritsky 2009; Olsen 
1999] |Holtzman et al.|[l999l |Dolphin||2000| |Smecker-Hane et al.| 
2002; Subramaniam 2004; Javiel et al. 2005 1. Furthermore, the star 



cluster population of the Large Magellanic Cloud has been inves- 
tigated, and ages and masses of a large fraction of the star clusters 
have been determined ijPietrzynski & Udalski 20 00] |Hunter et af] 



|2003| |de Grijs & Anders|2006 1. This enables us to study the star 
formation history of a galaxy with two independent methods. 

These introductory remarks outline the structure of this work, 
in summary: After first discussing the star cluster data set, we de- 
rive the star formation history of the Large Magellanic Cloud from 
the most-massive clusters (Sec |2.2} and from the total population 
(Sec. |2.3| l. Then we describe the results obtained from colour mag- 
nitude diagrams (Sec.[3j. We finish with a comparison of the results 
(Sec.|4]l and a summary (Sec. [5}. 



2 THE STAR FORMATION HISTORY OF THE LARGE 
MAGELLANIC CLOUD AS SEEN BY STAR CLUSTERS 

2.1 Data 

The ages and masses of the star clusters we use for the analysis 
are taken fro m|de Grijs & Anders| < [2006| >, which re-analysed the 
photometry of Hunter et al. O003| l, which itself is based on the ob- 
servations by Massey (2002f. The rectangles in Fig. [T] show the 
spatial coverage of the observed regions with star clusters. The 
star cluster ages and masses were derived by |de Grijs & Anders] 
(20061 from broad-band spectral energy distributions using their 
AnalySED tool which is based on the GALEV single stellar popu- 
lation models jKurth et al.| 1999, Anders & Fritze-v. Alvensleben| 
|2003| |Anders et aT]|2004^ The age uncertai nties for the 922 sta r 
clusters are in the range Alog,g(x/yr) ^ 0.35 (de Grijs et al. 2005 l. 
Due to discrete isochrones the age-mass diagram in Fig. |2 shows 
columns of star clusters of the same age. The lower mass limit for 
detection of clusters increases with increasing cluster age, as clus- 
ters fade due to stellar evolution and dynamical loss of stars, lead- 
ing to the wedge-like shape of the data in Fig. [2] 

The AnalySED tool provides "initial" masses of the star clus- 
ters which are corrected for mass loss due to stellar evolution. How- 
ever, the mass of a star cluster diminishes in time also because stars 
are lost in consequence of dynamical evolution. As we need true 
initial masses for the star clusters we correct for the dynamical 
evolution using the formulae of La mers et al.| ( |2005) . Given the 
dissolution time of a 10 4 M Q star cluster, f4, the initial star cluster 
mass of age t is given as 



Mi, 



M(t) 
M 



0.62 



0.62 x t 
O4/66O) 1 034 



1.61 



(1) 
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(This follows from combining eqq. 7 and 11 of Lamers ~et al.| 
(20051, with y= 0.62 and omitting the term for stellar evolution 
in eq. 11). For the Large Magellanic Cloud Boutloukos & Lamers 
1 2003} fou nd log 10 ?4 = 9.7 (using a smaller data set) and |de Grijs| 
|& Anders 1 2006 1 gave the slightly larger value of log 10 tt\ = 9.9 (8 
Gyr). |Parmentier & de Grijs](2008F carefully performed a reanaly- 
sis of the dissolution time and concluded that it is with the current 
data set only possible to constrain Ia, to be larger than 1 Gyr. There- 
fore we use two values for fy, 1 Gyr and 8 Gyr, to correct for dy- 
namical evolution. To visualise the difference between these values 
we show in Fig.[2]the not back-evolved masses as grey circles, and 
the "true" initial cluster masses as dots, using t\ = 8 Gyr in the top 
panel and ti\ = \ Gyr in the bottom panel. For us a larger ?4 seems 
to be more realistic, as the Small Magellanic Cloud has a similar 
value (log 1() f4 = 9.9, |Lamers et al.|2005] >, and more massive spiral 
galaxies with a deeper gravitational potential have smaller values. 
We will, however, discuss below the implications of both values 
when determining the star formation history from the most massive 
clusters. 

Further features in the age-mass diagram besides the typical 
wedge-like shape were pointed out by de Grijs & Anders (20061: 

(1) The large densities of clusters at log 10 x of 6.6 and 7.2: These 
are caused by the fitting procedure. There are no isochrones for 
clusters younger than 4 Myr (log 10 x = 6.6), and at logjQt = 7.2 
the isochrones are discrete due to rapid evolution. This does not 
have a large influence on the determined SFH. 

(2) The under-density of data points between ~ 3 Gyr and 13 Gyr 
(rs 9.5 ^ log 10 T ^ 10.1), which is the "well-known LMC cluster 
age-gap". 

(3) Overdensities at 7.8 sC log 10 x «C 8.0, 2.8 s: log 10 (Af/M@) sC 3.4 
and 8.2 ^ lgx ^ 8.4, all masses. These features could be caused by 
the last encounters between the Large and Small Magellanic Cloud, 
but this cannot be concluded with sufficient certainty because of the 
lack of better age resolution and lack of orbital information for the 
galaxies. 

It has also to be noted that this star cluster sample does not 
contain the 30Dor region, containing the young star cluster R136. 
It was classified as a newly formed star cluster ("NC") by |Bica| 
|et al. 1 1999 1 and so in a group of objects which were not selected 
by |Hunter*e t al. (2003). However, R135 is a massive star cluster 
having a mass of 5.5 x 10 4 M s ( Hunter et al.|l995 i, and is the 
most massive star cluster recently formed. The inclusion of this 
cluster is therefore crucial to the method used in the next Section. 

2.2 Star Formation History using the most massive star 
clusters 



In Maschberger & Kroupa ( 2007 ) we presented and tested a method 
to derive the star formation history of a galaxy using the most 
massive clusters. This method is based on the observation that the 
brightness of the brightest young cluster in a galaxy is correlated 
with the (present) star formation rate (Larsen 2002, Weidner et al. 
2004; Bastian 2008). This can be understood following the argu- 
ment of | Weidner et al.| (2004| >. Within a certain time span of the 
galaxy's lifetime, 8f, the amount of mass assembled in (long-lived) 
stellar clusters is proportional to the star formation rate, 

M c i U8te8 =ASFR8f (2) 

(A is the proportionality constant). This mass in clusters is related 
to a number of clusters that have formed, 



N, 



clusters 



-^clusters 

M ' 



(3) 



where a universal cluster mass function is assumed to calculate the 
average mass of a star cluster, M. Interpreting the star cluster mass 
function as a probability distribution, this then allows one to cal- 
culate the distribution of the most massive star cluster, M max , that 
would be expected for the given Adjusters ■ From this model follows 
a relation of the mass of the most massive star cluster with the star 
formation rate within St. This can be inverted to SFR = /(M max ) 
which can be used to determine the star formation rate over time, 
discretised by St. 

In general M max follows a probability distribution, related to 
the star cluster mass function, which has to be taken into account 
for the inversion (details of this can be found in Maschberger & 
Kroupa 2007 ). To minimise the number of assumptions, especially 
the exact form (pure power law or Schechter function as suggested 
by |Gieles et al.|20 06) and parameters of the cluster mass function, 
we use the relation of the mean mass of the most massive cluster 
and the star formation rate. The M max -SFR relation can be directly 
calibrated with the observed relation of the brightest young cluster 
and the star formation rate in a galaxy (assuming that the brightest 
cluster is also the most massive one of the most recent time inter- 
val, an assumption which is discussed in more detail below). This 
M max ~SFR relation is then applied to a mean mass of the observed 
most massive clusters over several St (choosing the number of used 
St such that during the whole time of averaging the star formation 
rate in the galaxy is not changing significantly). By using a moving 
averaging window (moved in steps of 8f ) the time resolution of the 
obtained star formation history can be increased. The length of the 
averaging window is essentially constrained by the age uncertain- 
ties of the star clusters, which are constant in logarithmic space, so 
that we keep the averaging window also constant in log I0 . 

By using M max and the empirical calibration we have avoided 
the need of the exact knowledge of the star cluster mass function. 
However, another crucial ingredient in this method is the formation 
epoch, 8f , which needs more explanation. In this context the often 
mentioned "size-of-sample" effect has to be discussed. The "size- 
of-sample" effect is simply the statistical increase of the mass of the 
most massive cluster with increasing sample size. With the general 
assumptions of an unchanging cluster mass function and constant 
cluster formation rate (number per time) a logarithmic age-mass di- 
agram has the characteristic upper envelope of an increasing mass 
with time. Equally-spaced time intervals in logarithmic space con- 
tain more physical time, thus more clusters are formed and sub- 
sequently the mass increases. However, this is not the full picture 
as the cluster (or star-) formation rate can change with time, lead- 
ing for example to the "age gap" in the Large Magellanic Cloud 
where barely clusters are found. A mathematically more correct de- 
scription would be a star cluster mass function depending on both 
mass and time, which is however not very practical. Here the for- 
mation epoch, 8? , comes into the play: this is the time by which the 
time evolution of a galaxy is discretised. With a reasonable choice 
of 8f the star formation rate in the galaxy can be assumed to stay 
constant, simplifying the statistical treatment. The increasing enve- 
lope in the log(age)-log(mass) diagram is preserved with using St 
(shown in fig. 3, top panel, of Maschberger & Kroupa 2007J. The 
difference to the established understanding of the "size-of-sample" 
effect is that one does not increase the size of a single sample, but 
one instead increases the number of samples. 

The question is now what a reasonable size for 8f is. Already 
mentioned was the need for the star formation rate to be constant 
over St, which gives an upper limit for 8? of ~ 100 Myr, the dy- 
namical time of a galaxy. In fact, the star formation rate should 
be constant over several St so that it can be averaged over sev- 
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Figure 3. Star formation history of the Large Magellanic Cloud derived 
from the most massive clusters. In the top panel we shows the age-mass 
diagram of the star clusters (with dynamically back-evolved masses, (4=8 
Gyr), including 30Dor as the star and highlighting the most massive clus- 
ters of each formation epoch with bigger circles. The dashed line is the 
fading limit. The lower panel shows the star formation history derived from 
the most massive clusters (with back-evolved masses and f4=8Gyr) as thick 
solid lines, splitting into two branches, the upper using the fading limit mass 
in gaps and the lower using M max = 0M . . The thick dashed branch at early 
ages follows with including 30Dor. The gray area is the uncertainty propa- 
gated from the uncertainty in the cluster masses. A horizontal line marks the 
average star formation rate and is enclosed by the la and 2c curves (solid 
and dashed) derived from the statistical spread of the M milx . 



eral M max . We follow here the choice of Weidner et al. ( 2004) and 
|Maschberger & K roupa (2007 ) of 10 Myr, for the practical reason 
that most clusters are of this age in the observational M max — SFR 
plot, as the luminosity of a star cluster peaks at about 10 Myr. Thus 
the brightest cluster is in many cases of this age and at the same 
time the most massive. With a different choice of 8t the bright- 
est cluster would have to be replaced by the actual most massive 
cluster in the normalisation, i.e. a true M max -SFR diagram would 
have to be observed. One possible interpretation of Sf = 10 Myr 
would be that it is the typical time-scale on which the inter-stellar 
medium rearranges itself into a coeval population of star clusters 
that are distributed according to the star-cluster initial mass func- 
tion (cf. [Weidner et al. 2004 1. Anyway, the comparison of the star 
cluster results with the CMD results will give an indirect check if 
our normalisation is correct. 

The method of going in 10-Myr-steps through age to estimate 
the star formation rate only works if there are star clusters in each 
formation epoch. However, an inspection of the age-mass diagram 
(Fig. [2} shows that for ages larger than 1 Gyr there are long inter- 
vals without clusters. For a formation epoch in which no cluster is 
detected a direct estimate of the star formation rate is not possible. 
But a non-detection does one allow to estimate limits for the star 
formation rate. If no cluster is detected, then none has formed mas- 
sive enough to be detected, which implies a low star formation rate. 



By using a star cluster evolution model the observational limiting 
magnitude can be translated into a mass, the fading limit. As no 
star cluster is detected above the fading limit, all cluster that might 
have been formed in this formation epoch must have had smaller 
masses. Thus the fading limit can be used as an upper limit for the 
most massive cluster in this epoch, and with this mass an upper 
limit for the star formation rate can be calculated. The lower limit 
for the star formation rate in an empty epoch is no star formation. 

This treatment of empty formation epoch gives us two results 
for the star formation history, an upper and a lower limit. When the 
upper and lower limit are identical (i.e. when only "full" forma- 
tion epochs are used), the most massive cluster method (using the 
M max -SFR relation) gives an estimate of the star formation history, 
otherwise the star formation rate can only be constrained by an up- 
per and lower limit. This is an implicit quality assessment of the 
method. Further uncertainties of results are caused by the proba- 
bilistic nature of M max , which mainly affects younger ages, and the 
uncertainties in the masses. We discuss these and their treatment 
after presenting the results for the Large Magellanic Cloud. 

Besides the availability of data there are two other sources of 
uncertainty for the most-massive cluster method, statistical scat- 
ter and the age/mass uncertainties. Due to the probabilistic nature 
of M max the stochastic scatter in the recovered star formation rate 
is very large for young ages, as averaging occurs only over a few 
8f, and decreases with increasing time. The amount of stochas- 
tic scatter has been determined from Monte-Carlo experiments by 
Maschberger & Kroupa (2007, their sec. 4.3, eqq. (16) and (17)). 
The age uncertainties are accommodated for by the averaging win- 
dow, which is kept constant in logarithmic time and has approxi- 
mately the size of the age uncertainties (0.5 dex). The mass uncer- 
tainty is propagated to an uncertainty in the star formation rate by 
using the M max values plus/minus their uncertainty in mass. 

In Fig. [3] we show in the upper panel the age-mass diagram of 
the star clusters, where the clusters identified as M max are the large 
dots. R136 is shown as an open circle, as it is not contained in the|de] 
Grijs & Anders ( 2006) sample. The dashed line is the fading limit, 
the mass that a cluster with the lowest observed brightness would 
have (calculated with the GALEV models). The upper and lower 
limit for the star formation history are shown as bold lines, which 
have the same values up to ~ 1 Gyr (calculated with = 8 Gyr). 
For the youngest ages the dashed branch follows by including R136 
in the star cluster sample. The uncertainty in the star formation rate 
introduced by the uncertainties in the cluster masses is visualised 
as a grey region. To assess the significance in variations of the star 
formation rate we show a constant star formation rate (the thick 
solid line at ~ O.lM /yr) which is embraced by the statistical la 
and 2o" scatter (thin solid lines and dashed lines, see Maschberger 
|& Kroupa|2007l ). 

Generally the obtained star formation history follows the dis- 
tribution of the star clusters for about one Gyr, when the number of 
clusters starts thinning out. Both the upper and lower limit agree for 
that period, so that the result of the method should be an estimate 
of the star formation history until that age. The peaks in the star 
formation rate are somewhat displaced when compared to the loci 
of the massive clusters which is caused by the time averaging. For 
ages younger than 100 Myr it is not possible to establish whether 
the variations in the derived star formation rate are caused by varia- 
tions in the actual star formation rate of the Large Magellanic Cloud 
because of the large stochastical scatter. In the age range from 100 
Myr to 1 Gyr the derived star formation rate deviates ~ 2a from a 
constant star formation rate, which should be caused by a decrease 
in the actual star formation rate. 
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Figure 4. Influence of t4 on the star formation history of the Large Mag- 
ellanic Cloud using star clusters. For the solid line ?4 = 8 Gyr and for the 
dashed line ?4 = 1 Gyr was used. 



We now turn to discuss the effects of the two values for t^. Fig- 
ure [4] shows this, with the star formation history using t\ = 8 Gyr 
as the solid line and t$ = l Gyr as the dashed line. Our results use 
only massive clusters, which are significantly affected by dynami- 
cal evolution only after a long time. Therefore the two solutions for 
the star formation history differ only at large ages. Essentially, the 
shorter ?4 implies a stronger dynamical evolution of the clusters, 
which consequently had larger initial masses, leading to a higher 
derived star formation rate. The differences in the star formation 
histories for ages younger than ~ 1 Gyr are only small. For larger 
ages the star formation rates are by a factor of ~ 10 larger for the 
smaller f 4 . 



2.3 Star Formation History using the total mass in star 
clusters 

As the Large Magellanic Cloud is very near to the Milky Way not 
only high-mass but also intermediate-mass clusters (with masses 
down to a few thousand Mq) can be detected over an extended 
time span. This allows us to use not only the most massive clusters 
to derive star formation rates, but also the whole cluster population. 
The amount of star formation in star clusters in a given time interval 
is simply the ratio of the total mass of star clusters and the length 
of the interval, 



SFR(?,M C 



1 I 
At *-* 

Mj>M c 



(4) 



where Mf is the mass of the i — th cluster with its age t(Mj) and M c 
is the completeness mass which follows from the detection limit. 

We derive the history of star formation in star clusters by mov- 
ing a time interval of constant logarithmic size (0.35 dex) in steps 
of 1 Myr until its boundary reaches an age of 400 Myr. For older 
ages the number of star clusters is too small to reach reasonable 
results. We use dynamically back-evolved cluster masses, M,-, with 
f 4 =8Gyr. 

Because the observations do not reach down to the lowest 
masses which star clusters can have, SFR gives only a fraction of 
the total star formation rate in star clusters. Thus the mass in star 
clusters has to be extrapolated to a total mass of stars in star clus- 
ters by assuming a star cluster mass function, £ oc M~$, a power 
law parametrised by an exponent (p) and a lower and upper mass 
limit (Mi and M u ). ^ is here normalised as a probability density, i.e. 
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Figure 5. Recent history of star formation in star clusters, derived by 
adding up all cluster masses above a completeness mass, M c and re- 
normalised as described in the text (cluster masses were dynamically back- 
evolved using ?4=8 Gyr). R136, shown as a star, is not included in the analy- 
sis. Various values for M c have been chosen to demonstrate the dependence 
of the result on the completeness. The top panel shows the location of the 
constant (4000, 7000 and 10000 M a ) and time-variable M c (407, 646 and 
1023 Mr. at t = 0) in the cluster age-mass diagram. For the constant M c the 
star formation histories are shown in the middle panel. The bottom panel 
contains the star formation histories for time variable M c , also showing the 
results from the middle panel as grey lines for comparison. 



fraction of star clusters, 



M ohs 
as 

a(M c ) 



The correct star formation rate in clusters is then 



SFR(r,M c 



a(M c )SFR(tM c 



(5) 



(6) 



(7) 



1. The normalisation factor follows from the observed 



The derived star formation history depends on the chosen complete- 
ness mass and parameters of the star cluster mass function. 

In order to explore the robustness of the obtained results we 
turn first to the completeness mass, as it seems not to be too well 
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constrained in our data set (see e.g. the discussion in Par mentierj 
|& de Grijs|2008| and paschberger & Kroupa|2009) To circumvent 
this problem we choose various possibilities for M c , shown as lines 
in the top panel of Fig. [5] the age-mass diagram. For a minimal de- 
pendence on the parameters of the cluster mass function we choose 
M c constant in time, with different values. Therewith one prohibits 
systematic effects in the shape of the star formation history caused 
by wrong parameters for the cluster mass function as the normalisa- 
tion factor is constant for all age bins. The results are shown in the 
middle panel of Fig. [5] for a better comparison re-scaled such that 
the star formation histories lie near together (actually using in all 
cases (3 = 2, M\ = 100 M Q and M u = 10 6 M ). The general struc- 
ture is the same for all star formation histories, except for the very 
youngest ages because of a lack of massive clusters. Further, small 
variations appear around ~ 15 Myr, ~ 60-90 Myr and after ~ 300 
Myr. The first two small variations are caused by an insufficiently 
small number of clusters more massive than M c . The discrepancy 
after 300 Myr is almost certainly caused by a too low mass for the 
lowest M c , leading to an incomplete data set at these ages. 



For an optimal use of the available data we choose a time- 
variable M c , running parallel to the lower envelope of the star clus- 
ters in Fig. [5] top panel. Now each age bin has an individual normal- 
isation constant, potentially introducing time-dependent systemat- 
ics. The lower panel of Fig. [5] shows the obtained star formation 
histories with the results for the constant M c (i.e. the results of the 
middle panel) plotted in grey for comparison. The overall structure 
of the star formation history is the same as for constant M c , with 
the exception that the peak at 90-150-Myr is more like a plateau. 
The small variations at ~ 15 Myr, ~ 60-90 Myr disappear with the 
larger number of clusters used, but the feature at ~ 300 Myr is still 
present. The agreement for the different choices of M c is better than 
for constant M c , which is rather surprising as the influence of the 
star cluster mass function is changing over time. This indicates that 
the results are robust and no systematical effects are introduced by 
the time-variable M c . 



For the correction of SFR we used (3 = 2, M\ = 100 M and 
M u = 10 6 M©. These values, especially P, are chosen such that the 
different M c all lead to the same result. As there are various values 
reported for (3 in the literature (e.g. Maschberger & Kroupa 2009 
|Gieles|[2009| and references therein and Weidner et al. 2004) we 
show in the top panel of Fig. [6] star formation histories corrected 
with different values of P (1.6, 2.0 and 2.4). M c is variable in time, 
starting with 260 Mq, (which is the second to lowest line in the top 
panel of Fig. [5}. For larger p the fraction of star clusters below 
M c increases, so that the star formation histories start at higher star 
formation rates. The increase of the star formation rates for differ- 
ent p is also time-dependent for time-variable M c , so that for older 
ages the amount by which the star formation rates are corrected in- 
creases. This leads to the growing difference between the curves in 
Fig. [6] top panel. The overall structure, however, remains the same 
within our range of P, and no additional features are introduced. 



The lower limit of the cluster mass function only has minor in- 
fluence on the absolute level of star formation histories, as evident 
in the lower panel of Fig. [6] Here we varied M\, using M\ = 5 Mr.,, 
50 M and 500 M (P = 2.0 and M u = 10 6 M ). The star forma- 
tion rates are a factor of 1.6 higher for M\ = 5 M :, compared to 
their values for M\ = 500 M Q . 
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Figure 6. Influence of the parameters of the star cluster mass function on 
the normalisation of the star formation history, derived from the complete 
cluster population (dynamically back-evolved cluster masses with ?4 = 8 
Gyr and time-variable completeness mass M c , the middle in Fig. [5] top 
panel). The top panel shows the changes caused by different exponents p" 
of the cluster mass function, especially on the older-age history. The choice 
of the lower limit M/ (bottom panel) has only minor influence on the abso- 
lute value of the star formation rates. x 
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Figure 7. Star formation history derived by Harris & Zaritsky 1 2009) using 
the colour-magnitude diagram method. The solid line is the SFH for the 
regions with in the Massey fields only (open circles in Fig. [TJ, with the 
uncertainty given by the grey region. The dotted line above is for the whole 
of the Large Magellanic Cloud. 



3 THE STAR FORMATION HISTORY DERIVED FROM 
COLOUR-MAGNITUDE DIAGRAMS 

For the comparison of our star cluster results with the results utilis- 
ing colour-magnitude diagrams we use the work of Harris & Zarit- 
|sky|p0 09). They presented the star formation history derived with 
the StarFISH software ( |Harris & Zaritsky 2001 ) in a grid of fields 
covering the whole Large Magellanic Cloud (the coverage is shown 
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Figure 8. Comparison of the results for the star formation history in the 
Large Magellanic Cloud. The dashed line within the grey shaded region 
(its uncertainty) is derived from colour magnitude diagrams for the parts 
lying in the Massey 1 2002 1 fields (as in Fig. ^] Sec. [3J. The thick solid 
lines are the solution using the most massive clusters, with the upper and 
lower limit for older ages, assuming no cluster formation or the maximum 
non detectable cluster mass in empty age bins (Fig. [3] bottom panel, Sec. 
|2.2) . The lowest dotted line is derived from the total star cluster mass (time 
variable M c (middle choice from top panel of Fig. [5] M; = 100 M s , fj = 
2, for details see Sec. \23) . Cluster masses have been dynamically back- 
evolved using t4 8 Gyr. 

in Fig.[TJ. Their photometric catalogue contains 24 million objects, 
so that each of the individually analysed fields contains some 10 4 
stars. For the synthetic colour magnitude diagrams the isochrones 
of the Padova group were used (Girardi et al. 2002;). The temporal 
resolution of the star formation history is given by age bins of 0.3 
dex size for ages younger than 100 Myr and bins of 0.2 dex for 
older ages. Solutions for the star formation history were obtained 
for four metallicities, Z = 0.001, Z = 0.0025 (interpolated and used 
only for ages larger than 100 Myr), Z = 0.004 and Z = 0.008. 

As their photometry does not reach the main sequence turnoff 
point for the old population the extraction of the early star for- 
mation history was difficult. Therefore |Harris & Zantsky| p009) 
restricted StarFISH to fit only a single age bin covering all ages 
older than 4 Gyr in the bar region. Within the bar region they used 
the typical star formation history from solutions for the star forma- 
tion history derived using HST data jQlsen|1999||Holtzman et~al"| 
1 1 999| and[Smecker-Hane et al. 2002] which widely agree with each 
other). 

Figure [7J shows the star formation history, derived by |Harris| 
|& Zar itsky ( 2009]). The star formation history of the fields within 
the regions of Massey (2002 ) (open circles in Fig. [TJ in which star 
clusters have been observed is shown as the solid line, with the 
grey area being its uncertainty. This partial star formation history 
follows the total star formation history (for the whole area of the 
Large Magellanic Cloud, dotted line) at about half the star forma- 
tion rate. 



4 COMPARISON OF THE METHODS 

In Figure [5] we summarise the solutions for the star formation his- 
tory of the Large Magellanic Cloud derived from colour-magnitude 
diagrams, the most massive star cluster or the total star cluster 
mass. As before in Fig. [TJ the dashed line within the grey shaded 
area gives the star formation history with uncertainty derived from 
colour-magnitude diagrams within the Massey (2002) fields. The 



thick solid lines give the upper and lower limits for the star for- 
mation history derived from the most massive star clusters (Sec. 
\2.2\ and the dotted line gives the amount of star formation in star 
clusters (Sec. [23}. 

The comparison of the M max estimate for the star formation 
history with the result from the CMD shows good agreement for 
ages up to ~ 1 Gyr, choosing the M max solution which includes 
R136. The absolute value for the star formation rate, derived from 
the most massive clusters, is at the level of the CMD star formation 
rate within the Massey (2002) fields. However, as the most mas- 
sive cluster-method is intended to give the star formation rate of an 
entire galaxy, and the spatial coverage of the star clusters contains 
most of the area recently active in star formation, it is perhaps more 
appropriate to compare to the galaxy-wide star formation rate. In 
this case the star formation rate would be underestimated by a fac- 
tor of 2 (compared to the dotted line in Fig.|7j and the normalisation 
in eq.|2]would need adjusting. As long as the spatial coverage is in- 
complete it is unfortunately impossible to disentangle inappropriate 
normalisation and effects of spatial incompleteness. 

Compared to the CMD solution the most massive cluster solu- 
tion shows an offset of the peaks in star formation. This is perhaps 
accounted for by the moving averaging window, although different 
isochrone sets could also account for this. If the ages of the star 
clusters could be more accurately determined the moving window 
could be reduced, by at the same time including not only the most 
massive star cluster but also the second, third etc. most massive and 
thus keeping the sample of data points used large enough for the 
statistically necessary averaging. An investigation in this direction 
is beyond the scope of this work. 

We turn now to a the comparison of the star formation rate 
derived from the total mass with the CMD and M tot solution. For 
the very youngest ages (< 10-20 Myr) the M tot star formation his- 
tory shows an increase of the star formation rate with age, in con- 
trast to the decrease of the CMD solution. This is because very 
young clusters are in a class of objects not selected by Hun ter et al.| 
( 2003 ) and thus the used data set is incomplete for very young ages. 
Within the age range where we can assume that the M tot solution 
is based on a complete data set (both in object selection and lower 
mass completeness), from ~ 10-20 Myr to ~ 200 Myr, the M to t 
solution shows the same structure as the CMD and M max solution, 
but at a lower star formation rate. The fraction of star formation 
in star clusters, i.e. the ratio between the CMD curve and the M to t 
curve, appears for the whole age range to be at a 10-20% level. 
This is caused by either the formation of only a fraction of stars in 
star clusters, or, when assuming that all stars form in custers, by 
the dissolution of star clusters caused by the presumably violent 
transition from their gas-embedded to the gas-free state. The sec- 
ond explanation could in principle be detected in the comparison, 
it should lead to a M tot star formation rate that is identical to the 
CMD star formation rate at very young ages (the embedded clus- 
ter phase). However, as the data set is not complete at the youngest 
ages a distinction between the alternative explanations is not possi- 
ble. 

For ages older than 1 Gyr the number of detected star clus- 
ters is very small, and their distribution shows many gaps. Thus no 
large amount of information about the star formation history can 
be extracted from the star clusters. The M max method gives an up- 
per and lower limit for the star formation history, either assuming 
no star cluster formation in the gaps (lower bold line in Fig. [8} or 
assuming the mass of the detection limit as the upper limit for the 
mass of a cluster that could have formed in an gap (upper line). As 
the mass associated with the detection limit increases with increas- 
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ing age, the upper limit for the star formation history shows also an 
increase of the star formation rate with age. The uncertainties intro- 
duced in the M max solution by stochastical scatter and mass errors 
are small compared to the effect of missing data. As there are some 
clusters with very high initial masses and ages around 10 Gyr one 
would conclude that there was a very high star formation activity, 
which is visible in the lower limit provided by the M max solution. 

Given the good agreement of the M max and CMD method at 
younger ages one would expect that the CMD solution falls in the 
region between the two M max limits. However, the CMD solution 
shows at an age of 2 Gyr a peak in the star formation rate which is 
nearly a factor of 10 higher than the upper limit for the star forma- 
tion rate following form M max . Although there are some clusters at 
2 Gyr with ~ 10 5 M (which lead to some kind of "peak" there in 
the lower limit of the M max method), many more would be needed 
to produce a signal. It also seems to be odd that there are essen- 
tially no clusters with masses between 10 4 M (the observational 
completeness limit, cf. Fig. [3} and 10 5 M Q (the mass of the ob- 
served clusters), a mass range which should be populated assuming 
a normal cluster mass function. 

If star formation is coupled to star cluster formation, as as- 
sumed for the M max method, and as shown by the agreement for 
ages < 1 Gyr, there should actually be also a large number of star 
clusters with ages of ~ 2 Gyr. There are perhaps several explana- 
tions for these missing clusters: they are located outside the ob- 
served fields, or a very different mode of star formation was active, 
where no massive clusters are produced, or a cluster destruction 
mechanism which only affected the missing clusters (e.g. interac- 
tions between the Large and Small Magellanic Cloud as suggested 
by |Bekki et al.|2004l >. 



5 SUMMARY AND CONCLUSIONS 

We compared the results for the star formation history of the Large 
Magellanic Cloud derived from its star cluster population, either 
following the method of Maschberger & Kroupa (2007), using the 
most massive clusters only, or using the total mass of the whole 
mass range of clusters (but then only for the most recent 400 Myr), 
with the results derived from colour-magnitude diagrams by Harris 
|& Zaritsky|p009) . 

We found that the results using the most massive clusters, both 
the absolute value for the star formation rate and the structure of 
the star formation history, agree well for the first Gyr with the star 
formation history derived from a CMD. For ages larger than 1 Gyr 
only a small number of clusters is detected, so that only a lower 
and upper limit of the star formation history can be given with the 
M max method. The CMD solution does not fall within these limits 
between 1 and 3 Gyr, but shows a peak with a higher star formation 
rate. This implies that the number of detected star clusters is too 
small compared to the expectations from the star formation rate 
following from CMDs. One possibility to resolve this discrepancy 
would be that additional clusters are contained in the area which is 
not observationally covered. 

Furthermore we derived the star formation history using all 
available star clusters, which, however, is only feasible for the most 
recent 20-400 Myr. The shape if this star formation history agrees 
with the CMD and M max results, albeit with absolute values for 
the star formation rate a factor ~ 10 lower. This implies that the 
fraction of star formation in (presumably bound, open) star clusters 
after gas expulsion is at a 10-20% level. Alternatively, this means 



(assuming that all stars form in a clustered way), that star clusters 
have an infant mortality of 80-90%. 

Our results show that star clusters are a powerful means to 
investigate the star formation history of a galaxy and invite further 
investigation in that direction. 
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